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ly-^ . In the framework of a recently proposed topological approach to phase transitions, some sufficient 

f"**^ ■ conditions ensuring the presence of the spontaneous breaking of a Z2 symmetry and of a symmetry- 

breaking phase transition are introduced and discussed. A very simple model, which we refer to 
, as the hypercubic model, is introduced and solved. The main purpose of this model is that of 

illustrating the content of the sufficient conditions, but it is interesting also in itself due to its 
simplicity. Then some mean-field models already known in the literature are discussed in the light 
of the sufficient conditions introduced here. 
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I. INTRODUCTION 



Phase transitions are very common in nature. They are sudden changes of the macroscopic behavior of a natural 
^ ■ system composed by many interacting parts occurring while an external parameter is smoothly varied. Phase transi- 
1 t \ tions are an example of emergent behavior, i.e., of collective properties having no direct counterpart in the dynamics 
or structure of individual atoms The successful description of phase transitions starting from the properties of 
the microscopic interactions (3^ between the components of the system is one of the major achievements of equi- 
librium statistical mechanics. From a statistical-mechanical point of view, in the canonical ensemble, describing a 
system at constant temperature T, a phase transition occurs at special values of the temperature called transition 
points, where thermodynamic quantities like pressure, magnetization, or heat capacity, are non-analytic functions of 
£j ' T; these points are the boundaries between different phases of the system. Starting from the celebrated solution of 
the two-dimensional fsing model by Onsager [2j, these singularities have been indeed found in many models, and later 
developments like the renormalization group theory J| have considerably deepened our knowledge of the properties 
of the transition points, at least in the case of continuous transitions, or critical phenomena. 
. Yet, the situation is not completely satisfactory. First, in the canonical ensemble these singularities occur only 
_, ' in the rather artificial case of infinite systems: following an early suggestion by Kramers [j], Lee and Yang [5j 
showed that the thermodynamic limit N — > 00 (N is the number of degrees of freedom, and the limit is taken at 
fixed density) must be invoked to explain the existence of true singularities in the canonical partition function Z(T) 
and then in the thermodynamic functions defined as derivatives of Z(T). Since in the last decades many examples 
of transitional phenomena in systems far from the thermodynamic limit have been found (e.g., in nuclei, atomic 
clusters, biopolymers), a description of phase transitions valid also for finite systems would be desirable. Second, 
while necessary conditions for the presence of a phase transition can be found (one example is the above-mentioned 
need of the thermodynamic limit in the canonical ensemble), nothing general is known about sufficient conditions: no 
general procedure is at hand to tell if a system where a phase transition is not ruled out from the beginning does have 
or not such a transition without computing Z: only for some particular systems or class of systems one can devise ad 
hoc procedures. This might indicate that our deep understanding of this phenomenon is still incomplete. 

These considerations motivate a study of the deep nature of phase transitions which may also be based on alternative 
approaches. One of such approaches, proposed in Ref. @ and developed later 0, is based on simple concepts and 
tools drawn from differential geometry and topology. The main issue of this new approach is a topological hypothesis, 
whose content is that at their deepest level phase transitions are due to a topology change of suitable submanifolds of 
configuration space, those where the system "lives" as the number of its de gree s of freedom becomes very large. This 
idea has been discussed and tested in many recent papers [cj Wi fTol Ull IT^ . Tl3L IbH ITU ITtj Il7l ll8L . Moreover, the 
topological hypothesis has been given a rigorous background by a theorem 20] which states that, at least for systems 
with short-ranged interactions and confining potentials, topology changes in configuration space submanifolds are a 
necessary condition for a phase transition. However, the converse is not true (there are topology changes which are 
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not connected with a phase transition |l2j\ and no sufficient topological conditions have been obtained yet. The 
problem of finding the sufficient conditions for a phase transition remains one of the fundamental open problems in 
this field: an answer to this problem would also make the topological approach an ideal candidate to define phase 
transitions in finite systems, for topology changes in the relevant submanifolds of configuration space do occur in finite 
systems, so that a criterion to discriminate the "good" ones would make them a natural extension of the concept of 
a phase transition to finite N case. 

The present paper aims at contributing to the search for topological sufficient conditions by considering, instead 
of the general problem of a generic phase transition, the more particular - but still very general and important from 
a physical point of view - case of the spontaneous breaking of a discrete symmetry. Although a phase transition is 
a far more general phenomenon, which is linked in general with the breaking of ergodicityj^, l2l| and may or may 
not be accompanied by the spontaneous breaking of a symmetry, many interesting phase transitions do occur in 
nature via the breaking of a symmetry. One of the most familiar cases is ferromagnetism: in the ordered phase the 
(continuous) rotational 0(3) symmetry of Heisenberg magnets is spontaneously broken. As to discrete symmetries, 
the paradigmatic example is the Ising model on a lattice, or, if one wants to consider continuous variables, the lattice 
tp 4 model, where if d > 2 a ferromagnetic transition exists and is accompanied by the breaking of the global Z2 
symmetry of the Hamiltonian, i.e., the symmetry under the simultaneous reversal of all the variables (the two- valued 
spins Si for the Ising model or the continuous real variables for the lattice ip model). 

In the following we will consider only the case of a Z2 symmetry, even if we believe that it should be possible to 
extend our results to general discrete symmetries. We shall also restrict ourselves to systems described by continuous 
variables, because the topological approach can be defined only for these systems, even if in some cases we may refer to 
Ising-like discrete spin systems for illustrative purposes. As we will show in Sec. lllll under fairly general assumptions 
it is possible to state a sufficient condition for the presence of a Z2-symmetry breaking phase transition essentially in 
terms of the topology of the equipotential hypersurfaces in configuration space, provided some additional conditions 
on the behavior with TV are satisfied. Before stating and discussing this results, in Sec.|n]we will discuss at a general 
level the problem of the breaking of a Z2 symmetry and the basis of the topological approach. Then, after having 
stated the above mentioned result, in Sec. II VI we will illustrate it introducing a simple abstract model, and in Sec. IVl 
we will discuss them in the light of some physical models already studied in the literature. We will end with some 
concluding remarks in Sec. IVII 



II. GENERAL PICTURE 



A. Phase transitions with Z2 symmetry breaking 

At a qualitative level, the physical mechanism underlying the spontaneous breaking of a Z2 symmetry is quite well 
understood |22| . To begin with, let us consider a simple example with a single degree of freedom: a particle of unit 
mass in a double- well potential V(q), such that V(q) — V(— q), at a fixed temperature T. The dynamics of the particle 
will be described by a Langevin equation 

. dV 

q = -jq-~r + v(t) , (i 

dq 

where r](t) is a ^-correlated white noise whose amplitude is related to the friction coefficient 7 and to the temperature 
by the fluctuation-dissipation theorem, 

(n(t) V (t + t)) = 2 7 T<5(t) . (2) 

The "magnetization", i.e., the order parameter of the system, whose nonzero value signals the breaking of the Z2 
symmetry, is the time average of the position q. As long as the temperature is low, the dynamics of the particle is 
essentially an activated process, with two widely separated time scales: one small scale in which the particle oscillates 
in one of the two wells, and a large time scale r in which one observes jumps from one well to the other. On time 
scales t <C r, the symmetry appears to be broken, because the particle is confined in one of the two wells and the 
finite-time order parameter is nonzero: 



t 

q(t') dt' ^ t < t . (3) 



n 



According to Kramers' theory [2J 



r oc exp \^r) . (4) 
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where we have set ks — 1 as we shall always do from now on, AE is the height of the energy barrier the particle has 
to overcome in order to jump from the minimum of one well to the other, and is a finite quantity if V is the potential 
energy of a single particle, so that, even if the timescale r increases exponentially while decreasing the temperature, 
the order parameter (q) vanishes for any finite temperature: 

1 r* 

(q) = lim - / q(t') dt' = VT > , (5) 

t-«X> t Jq 

and symmetry breaking is possible only at T — 0. Note that increasing AE does not change the situation, unless one 
takes the limit AE — > oo where the symmetry is broken for any value of T. In any case, no phase transition between 
a symmetric and a broken-symmetry phase is allowed in this single-particle system. 

Nonetheless, we are interested in many-particle systems: it is just the number of degrees of freedom, N, which plays 
a fundamental role to make a symmetry breaking possible 22] . The potential energy V(qi, . . . , gjv) is now a function 
of N variables, still Z 2 -symmetric, i.e., V(q\, . . . , gjv) = V{— q%, . . . , ~qN)- The potential energy necessarily has two 
equivalent absolute minima related by the symmetry, but multidimensionality of configuration space means that there 
could be many possible routes to go from one minimum to the other. If now we denote by AE the minimum barrier 
to jump in order to connect the two minima, for low enough T Eq. |0J still holds 33] . This means that, if AE grows 
with N, in the thermodynamic limit the equilibration time scale r becomes infinite and the system is trapped in one 
of the two wells even for infinite times: the order parameter ^(X^i 1i) * s now finite, 

1 / N \ 1 r* N 

n (E *) - & A, m I § dt ' * VT < Tc ' (6) 

and the symmetry is broken for finite temperatures below a critical temperature T c . Note that in Eq. © the two 
limits do not commute, i.e., we must first let the system go to the thermodynamic limit and then take the infinite-time 
averages. Doing the other way round we would not get any transition, and no symmetry breaking would be present 

M 

The reason why this happens only for T < T c and not for any T is due on the one hand to the fact that the 
separation of timescales, and thus the activated process picture, holds only for sufficiently low temperatures, and on 
the other hand to that the above discussion is oversimplified: we have neglected the role of entropy, and the high 
dimensionality of the configuration space ensures that for sufficiently high T"s the entropy will always disorder the 
system. The above argument could be made more stringent using free energy barriers instead of energy barriers. 
However, our purpose was only to show, using an intuitive dynamical argument, how in a many-particle system 
symmetry breaking can occur in the thermodynamic limit, not to prove that it does occur. We have thus seen that 
one of the basic ingredients for the possibility of symmetry breaking is, besides the thermodynamic limit, that the 
height of the minimum barrier to overcome must grow with N: here is where, among other factors, the dimensionality 
of the system comes in. To make a familiar example, in one-dimensional Ising systems AE is constant as N grows, 
because after having flipped one spin all the others can be flipped without any extra energy cost, while is proportional 
to \/N in a two-dimensional system, because to flip a whole region of spins the energy is paid at the perimeter of 
the region, whose length scales as L d_1 for a d-dimensional lattice of length L. We recognize here the Landau-Peierls 
argument |25j to prove the existence of a finite-temperature symmetry breaking in a two-dimensional Ising system, 
and indeed the physical content of this argument is the same of the "dynamical" argument above, as to the energy 
part. The Landau-Peierls argument is much more powerful because the entropic contribution, and then T c , can be 
estimated too, provided we can efficiently count the relevant configurations, which however limits its applicability to 
Ising-like systems. The dynamical argument is valid for general systems with continuous variables, but remains at a 
qualitative level and does not provide a clear sufficient condition for the presence of a phase transition. As we shall see 
in the following, it is possible to translate it into the topological language, which does allow to state such a sufficient 
condition. But before doing that let us review the basis of the topological approach. 

B. Basis of the topological approach 

As already mentioned in the Introduction, where also relevant references to original papers were given, the topolog- 
ical approach to phase transitions is based on the "topological hypothesis" that phase transitions are due to suitable 
topology changes in some submanifolds of configuration space defined by the potential energy function. Here we want 
to recall which is the basis of this approach. 

Let us consider a Hamiltonian system with N degrees of freedom and standard kinetic energy, described by the 
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Hamiltonian 

JV 

2 



1 N 

n = -J2Pi + v (n> ■■■><!*) > ( 7 ) 



where V(q) (from now on q = {q}jLi) is the potential energy and the g^'s and the p^s (i = 1, . . . , N) are, respectively, 
the canonical conjugate coordinates and momenta, and are continuous variables; q S M, where M is the iV-dimensional 
configuration space manifold, and V is bounded below on M. The conhgurational partition function of such a system 
can be written as (we omit the (h N N\) _1 normalization factor because it is irrelevant for our discussion) 

/>oo 

Z N (/3) = / d(Nv)e-P Nv u N (v) (8) 



where f3 — T 1 and u>n(v) is the density of states at (potential) energy per degree of freedom v, i.e., the Liouville 
measure of the isopotential hypersurface £„, 

u)N (v)=^ v )=J^-^, (9) 

where Y, v is defined as the iVv-level set of the potential energy V(q), 

E v ={qeM\V( qi ,...,q N ) = Nv} , (10) 

and dY, is the volume element on E„. We write u>n as 

un{v) = (a N (v)) N , (11) 

where oat(w) = exp [sjv(i')] and sn{v) is the configurational entropy per degree of freedom. Now we can write (without 
loss of generality, we assume that the absolute minimum of V is zero) 

poo 

Z N (fi) = N dve N{io&aNi - v) -P v , (12) 
Jo 

and then, as N gets very large, we can evaluate the integral over v in © by replacing it with the largest value of the 
integrand: 

Z N (/3) = const. e "[sup„(i°g«jv(«)-/J«] ( 13 ) 
where up to now we have only assumed extensivity, i.e., 

a(v) = lim a N (v) = [n(Z v )] 1/N (14) 

exists and is finite, which, due to the physical meaning of the density of states, amounts to requiring that the specific 
configurational entropy is well defined in the thermodynamic limit [34| ; then we can write 

Z N (P) — ATe-W) M ( E ) = ATe-W) / Jj= , (15) 
w - > °° " / s- tfJ) l|VV|| 

showing that the only relevant contribution to the partition function comes from a single isopotential hypersurface 
Su(/3)> where v((3) is the value of v which realizes the supremum in Eci. 1131 if we assume that the function (3v — log a(v) 
has a single minimum which does not coincide with the extrema of the interval of definition, v{0) is the solution of 
the saddle-point equation 

^^-loga(u)]=0, (16) 
and coincides with the expectation value of V/N, 
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This amounts to saying that as N — » oo the support of the equilibrium measure reduces to the equipotential hyper- 
surface S-(^). In other terms, when computing the canonical ensemble average of a function A(q), we can write as N 
gets very large 

^^-W)LM' (18> 

where j4| s is the restriction of the function A(q) to E^f/j) and 

m =Lwir (19) 

A major topology change in a family of manifolds, depending on a continuous parameter v, occurring at some v c may 
induce singularities in the ^-dependence of their volume, whence the basic idea of the topological hypothesis: the 
deep origin of a phase transition might be concealed in the way the configuration space is foliated in level sets of the 
potential energy function, for a sufficiently "strong" topology change in the S^'s or the M v 's at some v c might induce 
a phase transition at a temperature such that v c — -h>(V) because, as we have just seen, at very large N the measure 
concentrates on a single "slice" . 

An important remark is in order. All these results can be reformulated considering the submanifolds M v = {q £ 
M\V(q) < Nv} instead of the S„'s; the relation between these two families of submanifolds is £„ = dM v . The reason 
for the possibility of substituting the M v 's to the S^'s is in the fact that the Liouville measure of M v is the same 
as that of £„ when N — > oo, and that topology changes of the S^'s do occur simultaneoulsy with those in the M v 's 
apart from very particular cases. In applications, using the M v 's instead of the S^'s may be easier (for instance, in 
some cases Morse theory allows a direct calculation of the topology changes in the M v 's using the potential energy as 
a Morse function [HE!). 

As already noted in Sec. U the idea of the topological hypothesis has been discussed and tested in some particular 
models, and a theorem has been proven showing that - for a wide class of systems - topology changes in the X„'s 
are a necessary condition for a phase transition to occurr 20] . In the following Section we are going to show how a 
sufficient condition can be derived in the case of the spontaneous breaking of a Z2 symmetry. 



III. SUFFICIENT TOPOLOGICAL CONDITION FOR Z 2 SYMMETRY BREAKING 

Let us now consider a system of the class with a potential energy which is Z2-invariant. As we have shown 
in the previous Section, in the thermodynamic limit the canonical ensemble average of a function of the coordinates 
is given by Eq. I|18l) . Let us now consider, instead of a generic function A{q) 1 a function whose average is a order 
parameter for the Z2 symmetry breaking, i.e., 

^) = ^f>. (20) 

i=l 

By inserting Eq. (|20|l into Eq. (T5J, at a first sight we are led to conclude that no symmetry breaking is possible also in 
the thermodynamic limit, because all the hypersurfaces £„, being level sets of the potential energy function V, must 
have all the symmetry of the function V itself, and in particular the Z2 symmetry; hence the order parameter is 
zero for any value of T . 

This conclusion is wrong. The £„'s may be composed of two or more disjoint connected components, and although 
the whole £„ must respect the invariance, each single connected component need not to be Z 2 -invariant: each one may 
be the image of another one under the symmetry operation. If £„ is made up of disjoint connected components, then 
the definition H18|) of the ensemble average is perfectly legitimate, but cannot be consistent with the actual behaviour 
of the system, because the representative point of the system can explore only one of the connected components. This 
may be easily seen if we think of the dynamics. Saying that when N — ► 00 the support of the measure is ^v(p) is 
equivalent to say that the measure is formally the standard Boltzmann weight 

e„(g;/3) = ^e-™-<«« (21) 

with effective potential Vqo given by 



(22) 
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at t = 0, the system will be in one of the disjoint components, and it will remain there forever: it can never jump to 
another one because this would require to jump over an infinite energy barrier. Hence, if £„ = U U • • • U £™ 
with EjnEj =0Vo,6, the restricted ensemble average being equal to a time average will be given by one of the 
following 

(A) a (P) — » = -J— [ A ^ ad ^ , a=l,...,n, (23) 

v ' vm a^~ z a wM m ii w ii 

where the correct value of a will be specified by the initial conditions, and 

If neither £" nor A(q) are Z2-invariant, definitions (|I8|I and (|23[1 do not yield the same result; in particular, the 
average of X^iLi 3» according to Eq. 1)23(1 may give a nonzero result, so that symmetry breaking is possible. 

We note that the disjoint connected components £" of E„ play a role analogous to that of pure states in the standard 
approach to the formal treatment of symmetry breaking |22,[2g, while their union (i.e., the whole £„), plays the role 
of the mixed state. 

We can now state the following 

Theorem 1 (sufficient topological condition for Z2 symmetry breaking) Let us consider a system of the 
class with N degrees of freedom and a potential energy V bounded below which is ^-invariant. Let the en- 
tropy per degree of freedom be well defined in the thermodynamic limit, i.e., the function a(v) defined in Eq. \1$ exist 
and be continuous and piecewise differentiable. Let £„ be the family of equipotential hypersurfaces of the configuration 
space M defined as in Eq. Without loss of generality, let min(V) = 0. Let v" > v' > be two values of the 

potential energy per degree of freedom V/N such that £„ = U" =1 £° Vu € (v',v"), with SJflEj = Va, b, and such 
that Va 3 b 7^ a : Z(Y,%) = £^ where Z is the ^-symmetry map on M, Z(q) = — q. 

Then, in the thermodynamic limit the 1i symmetry is spontaneously broken for all the temperatures T £ (T',T"), 
where T" > T > 0. 

Proof. Thanks to the hypotheses on the function a(v), the statistical average of the potential energy per degree of 
freedom, v(T), is given by Eq. (|I7|I and is monotonically increasing as T is varied from to +00, because 41 < 
would imply a negative (configurational) heat capacity, which is forbidden in the canonical ensemble | 25| . Then there 
exists T',T" > 0, T = v- y {v'), T" = v-\v"), such that v(T) £ {v',v") VT £ (T',T"). For the sake of clarity let 
us assume that for v £ (1/, v") the equipotential hypersurface is made up only of two disjoint connected components, 
= £j" U £~ Vf < v', with £j" n £~ = and Z(Ti^) — £^"; the extension to a larger number of components is 
straightforward. According to Eq. (|23|l . for T £ (T',T") the order parameter is 

and since the integrand is odd under Z, Z + = Z~ - this follows from the definition in Eq. I|24|l - and iJ(£+) = £~ 
we have 

m+ = -m~ ^ , (26) 

because in order to have rrv^ = each of the £^'s should be symmetric around q ~ 0, but this is impossible if they 
are disconnected and Z2-invariant. Then the symmetry is spontaneously broken when T < T' . □ 

Corollary 1 Let us consider a system of the class with N degrees of freedom and a potential energy V bounded 
below which is ^-invariant, an let it have two degenerate distinct absolute minima q + = — q~ whose value is min(V) = 
0. Let the entropy per degree of freedom exist as in Theorem^ Let Let v' > be a value of the potential energy per 
degree of freedom V/N such that all the £„ 's are homeomorphic for v < v' . 

Then, in the thermodynamic limit the Z2 symmetry is spontaneously broken for all the temperatures smaller than a 
finite value T' . 

Proof. When v — ► 0, £„ = £+ U £,7, where £^ are homeomorphic to (N — l)-spheres and £+ n £,7 = 0. Since no 
topology change occurs until v' , for all the temperatures smaller than T' = u _1 (w') we are in the situation of Theorem 
^ whence the thesis. □ 
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We note that the crucial ingredient to obtain the sufficient condition of Theorem ^ is that the value of the potential 
energy below which the topology of the E„'s is such that they are made up of disjoint connected components must 
be proportional to N, and this agrees with the qualitative reasoning of Sec. Ill Al Being a sufficient condition, but 
not a necessary one, this obviously does not rule out the possibility of having a symmetry breaking with less strong 
assumptions. Moreover, we have proven that the symmetry is broken below T", but neither that it is restored above 
T", nor that there is a phase transition at T" or at any other T. To prove that, we have to strengthen the hypoteses, 
as in the following 

Theorem 2 (sufficient condition for Z 2 -symmetry-breaking phase transitions) In addition to the hypothe- 
ses of Theorem^ let v'" > v" be a value of the potential energy per degree of freedom V/N such that £„ is made 
of a single connected component for v > v'" . Let also the support of the canonical measure be the whole when 
v(T) > v'" also in the thermodynamic limit. 

Then, in the thermodynamic limit there exist finite temperatures T' < T" < T'" such that the Z 2 symmetry is 
spontaneously broken for all T S (T',T") and is restored for all T > T'" . There is also at least a phase transition 
(or more than one) at T c such that T" < T c < T'" . If v" = v'" = v c , then there is just one phase transition at 
T c = w _1 (u c ). 

Proof. Theorem n ensures that m 7^ for T G (T',T") — (v _1 (v'), tJ _1 (ti")). Then, reasoning as in the proof of 
Theorem [U there exists T'" > 0, T" 1 = W 1 ^'") such that v(T) > v'" VT > T'". Then, for T > V" the order 
parameter is 



m 



N^^/n^oo NZ / v , W 



because the E^^'s are all symmetric: the symmetry is restored. Hence, the order parameter m is nonzero when 
T 6 (T',T"), and is constant and equal to zero when T > T'" . Then the function m(T) has at least a non-analytic 
point for a temperature T" < T c < T'", thus there is at least one phase transition in the system. If v" = v'" = v c , 
then there is a transition occurring at T c — v~ l (v c ). □ 

We note that the hypothesis that the support of the measure remains the whole Iwt) when v(T) > v'" also in 
the thermodynamic limit is essential, because without this assumption we cannot prove that m = 0. If as N — > 00 
the support of the measure shrinks to a submanifold of XWy) , then the symmetry may remain broken even if S^(T) 
is Z 2 -symmetric. We believe that this case is not a purely academic one: it is probably what happens in at least 
one physically relevant example (the mean- field ip 4 model, see Sec. |V}. This assumption was not necessary at all in 
proving Theorem ^ because this may affect only the actual value of m, but not the fact that it is nonzero. 

These two theorems have rather strong assumptions, which may probably be weakened. Moreover, for a generic 
many-particle system it is not an easy task at all to characterize all the topology changes undergone by the XVs 
(see e.g. the discussion given in Refs. 11 mm), so that we are not claiming to have derived a practical all-purpose 
method to prove the existence of symmetry-breaking phase transitions in generic systems. Nonetheless, the theorems 
do allow one to make predictions, which are confirmed by the analysis of some known models, to be discussed in 
Sec. But before doing that, let us discuss an abstract toy model we introduce in the next Section, whose main 
purpose is to illustrate the theorems in a simple and clear way, but which may have also some interest on its own. 



IV. HYPERCUBIC MODEL 



We now introduce and solve an abstract model to enlighten, in a pedagogical way, the content of the theorems 
proven in the last Section. This model is rather abstract, and from a physical point of view it can be seen as a model 
of a particle bouncing in a potential in an ./V-dimensional space, but we build it starting directly with the cquipotcntial 
hypersurfaces £„• 

The simplest Z2-invariant potential is a double square well in one dimension, i.e., 

( if a < \q\ < b; 
V(q) = \ v if a > \q\ ; (28) 
[ +00 if |g| > 6, 

with < a < b. In this simple case the configuration space M is just the real line K. Our first toy model, which we 
will refer to as the hypercubic model, is nothing but a generalization to N dimensions of this double square well. The 
configuraton space is now WL N , and in M we consider two disjoint hypercubes A + and A~~ , symmetric under Z 2 , and 
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FIG. 1: Sketch of the hypercubes A ± and B for N = 2. 



third hypercube B, centered in the origin, such that A + ,A~ C B. Then we define 



if q £ A ± ; 
V(q) = { Nv c if q £ B\{A+ Ui"} 
hoo if q£ R N \B. 



(29) 



This potential is Z2-invariant by construction. The hypercubes and B are sketched in Fig. Q] in the case N = 2. 
The equipotential hypersurfaces £„ are then 



if v < 0; 

U A~ if u = 0; 

E t , = { if < v < v c ; 

B\{A+UA-} if v = v c ; 

if v > v r . 



(30) 



Then the topology of the S„'s is such that (~ stands for "is homeomorphic to") 



D J 



if 


v < 0; 


N + D N if 


u = 0; 


if 


< w < v c 


£ 2 W if 


v = v c ; 


if 


v > «c , 



(31) 



where D N is a disk in R N , with we denote a two-punctuated disk (a disk with two disjoint disks removed) and 
"+" stands for the disjoint union. The fact that apart from the values v = and v = v c the S K 's are empty sets is 
due to the very singular nature of the potential, which has only two possible values instead of a continuous interval: 
as v increases starting from values smaller than zero, the system potential energy actually "skips over" all the values 
but and v c , because these are the only allowed values of V/N. Nonetheless, we see that the E„'s of the hypercubic 
model undergo a topology change as v changes from v = to v = v c of the kind described in Theorem [21 in the case 
v' = v" = v c . However, the situation is much clearer if we consider the M v manifolds instead of the S^'s: we have 



. if v < 0; 

M v = I A+UA- if < v < v c 



B 



if v > v c , 



(32) 



and the the topology of the M v 's is 



M„ 



if v < 0; 
d n + d n |f q < v < Vc . 

D N if v > v c • 



(33) 
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Then the M^'s of the hypercubic model undergo a topology change as v changes from v — to v = v c precisely of 
the kind described in Theorem [3 in the case v' = v" = v c , and Theorem [21 (applied to the M„'s) states that the 
hypercubic model, in the thermodynamic limit N — > oo, undergoes a phase transition with Z2 symmetry breaking at 
a finite temperature T c such that jr{V)(T c ) = v c . Let us see it explicitly, solving the model. 
At any finite N, the configurational partition function of the hypercubic model is 

Z N (p)= f d N qe-^=f d N q+ f d N q + e-^f d N q (34) 

JR N J A+ J A- J B\{A+UA-} 

and denoting by o and b the length of the side of A^ and B, respectively, we obtain 

Z N {{3) = 2a N + (b N - 2a N )e- l3Nv < , (35) 

where b > 2a because A + , A~ C B. Thermodynamic functions can then be computed at any finite N and their limit 
as N — * 00 can be studied directly. We shall see that in the following, but before doing that we note that when N is 
large we can write 

Z N ((3) — » 2a N + b N e-? Nv ° = 2e An ° ga + e N ^ b -^ , (36) 

so that in the thermodynamic limit only the largest of the two exponentials contributes to Z, and there will be a 
critical value j3 c of the inverse temperature /?, given by the equation 

log a = log b - (3 c v c (37) 

whose solution is 

p c = -\og(-\ , (38) 

such that 

7 (R\ — > / 2eJV1 ° SO if ^>^5 

This means that the system feels an effective double-well potential, with potential energy zero and the two wells 
separated by an infinite barrier, for (3 > (3 C and an effective single-well, symmetric potential with potential energy v c 
when (3 < (3 C . Hence the symmetry is broken when f3 > f3 c . The value of the order parameter 



, 1=1 



where the average is the restricted one when the symmetry is broken, will be, in the thermodynamic limit, 

where go is the value of all the coordinates qi,...,qx of the center of the A + hypercube. The order parameter is 
plotted as a function of T = /3 _1 in Fig. El I n the- same limit, the average potential energy per degree of freedom will 
be 

^=-W ogZN{0)= {l (42) 

so that the phase transition is a discontinuous (first-order) one. 

Computing the thermodynamic functions at finite N one finds for the average potential energy per degree of freedom 

v (b N - 2n N )e~P Nv " 

<"» ft (43) 
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FIG. 2: Order parameter of the hypercubic model in the thermodynamic limit, as a function of T. Here v c = 1, a = 1, b = 2a, 
so that T c = (log2)~ 1 . Only the positive branch is plotted. 



A 
> 
v 




FIG. 3: Average potential energy of the hypercubic model as a function of T. The different smooth curves are the finite-iV 
result Q430 with TV = 10, 20 and 50, while the piecewise constant curve is the N — > oo limit l|42|l . Numerical values as in Fig. [2] 



and for the configurational specific heat 



c v (p;N) 



2Nf3 2 v 2 c a N (b N -2a N )e-P Nv c 
[2a N + {b N - 2a N )e-f )N ^f 



(44) 



These functions are plotted in Figs. 01 and 0| We see that at any finite N both (v) and cy are regular functions, which 
converge (non uniformly) to the limiting non-analytic functions already determined before. As to the specific heat, 
since (v) is a piecewise constant function in the thermodynamic limit, cy = everywhere but at T c where it has a 
5-like singularity. 
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FIG. 4: As in Fig. |3]for the specific heat cy- 

The hypercubic model is undoubtely very abstract, for it describes a particle in a TV-dimensional potential which is 
highly singular: yet it is extremely simple and neatly illustrates some consequences of Theorem|21 Moreover, it can be 
probably considered as one of the most elementary models exhibiting a phase transition, so that it has a pedagogical 
interest on its own. 

However, as we shall see in the next Section, there are models, already known in the literature, which can be 
analyzed in terms of the theorems of Sec. lIIII and whose potential energy function is indeed regular and describes the 
interaction among microscopic degrees of freedom, although in a not completely realistic way due to its mean-field 
character. 



V. PHYSICAL MODELS 

A. Mean-field spherical model 

A physical example which can be exactly analyzed in terms of Theorem|5|is the mean-field spherical model, recently 
studied in Ref. [hHI^ . which is the mean-field version of the model originally introduced by Kac and Berlin in 1952 
[28j . It is an Ising-like system with continuous variables: the potential energy is 

1 N 

V&) =~^Y J WPi , (45) 
where the variables <p; G M, i = 1, . . . , N are subject to the condition 

N 

5>? = iV (46) 

8=1 

which constraints the variables to live on the (N — l)-dimensional sphere of radius ^/N centered in the origin of M. N , 
whence the name "spherical model" . The Z 2 invariance of the potential is apparent from Eq. (|45|) . The potential 
energy is also bounded, so that the potential energy per degree of freedom is bounded too: v € [— |, 0]. 

The topolo gy o f the equipotential hypersurfaces £„ of the mean-field spherical model can be easily determined, and 
one finds 0,H3 that 

!0 if v < -|; 

S N-2 +g N-2 if _1 < v < q. (47) 
S N - 2 if V = . 
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We are in the situation described in Theorem [2 which then predicts that this model, in the canonical ensemble, has 
a phase transition with Z 2 symmetry breaking at a finite tem per ature T c such that (v)(T c ) = 0. This is indeed what 
happens, and the critical temperature turns out to be T c = 1 [Tj|. 

The mean- field spherical model is thus a nice illustration of the consequences of Theorem |3 



B. Mean- field tp 4 model 

An example where Theorem^holds is provided by another mean-field version of a continuous-spin Ising model, the 
mean-field <p 4 model 0, 0, Il9j . The interaction potential is 

T N N / 1 i \ 

i,j=l i=l ^ ' 

where J > is a coupling constant and (fi € R, i = 1, . . . , N. Again the Z 2 symmetry is apparent from Eq. 1)480: 
this potential is bounded below, but not above, and v G [w m in, +00), where u m in = — 4 . There are two equivalent 
distinct minima, ipi = ±V J + 1, i = 1, . . . , iV, then as w — ► v m in 

E„ ~S Ar - 1 +S JV_1 . (49) 

It has been numerically shown [l^ that for any N and any J > there are no topology changes [3^ in the E„'s as 
long as v < v' , where v' < is a finite value which grows when J grows. Then, this model fulfills the hypotheses 
of Theorem or more precisely those of Corollary ^ so that the Z 2 symmetry must be broken below a finite 
temperature T' = tJ _1 (w'). This is precisely what happens: the magnetization m(T), which can be exactly calculated 
in the thermodynamic limit due to the mean-field character of the model, is nonzero for T < T' . For instance, as 
J = i, v' ~ —0.35 and T' ~ 0.3, while the critical temperature below which the symmetry is broken is T c ~ 0.4; when 
J = 1, v' ~ -0.25 and T' ~ 0.9, while T c ~ iJOl- 

However, as w > 0, it has also been shown |!7lll8lll9l | that 

E„ ~ S^- 1 , (50) 

but at least for sufficiently large values of J the Z2 symmetry remains broken also for temperatures T > T" = TJ _1 (0). 
Then, although this model is a clear illustration of Corollary ^ it does not behave as predicted by Theorem [2 at 
variance with the mean-field spherical model. This means that at least one of the hypotheses of Theorem El does 
not hold for the mean-field <p 4 model. As we will discuss elsewhere j^, in this case it is not true that the support 
of the equilibrium measure remains the whole E- even in the thermodynamic limit, due to the additional constraint 
provided by the fact that the function whose average is the order parameter enters the potential [37| 



VI. CONCLUDING REMARKS 



We have shown that the topological approach to phase transitions allows one to translate an intuitive, qualitative 
picture of the origin of discrete-symmetry-breaking phase transitions into a sufficient condition for such phenomenon 
to occur. This is a first step towards obtaining more general sufficient conditions for phase transitions, thus towards 
filling a gap in our present understanding of these ubiquitous and fascinating phenomena. The topological sufficient 
conditions have been derived here in the case of Z2 symmetry, but the proofs of Theorems ^ and [3 are easily adaptable 
to more general discrete symmetries: we chose to restrict ourselves to Z2 for the sake of simplicity, and because of 
the great importance of this particular symmetry in the development of our understanding of the physics of phase 
transitions. We note that our sufficient conditions are not exclusively topological in nature, i.e., they cannot be 
formulated in terms of the topological properties of the EJ,s alone: for example, we have also to introduce some 
hypotheses on the behavior with N of the values of the energy at which some topology changes must happen. 
Nonetheless, these additional requirements appear absolutely natural, and do not change the fact that the basic 
ingredient of the two theorems proven here comes from topology. 

A simple model introduced here, the hypercubic model, and some mean-field models already known in the literature 
have been discussed in the light of these two theorems. It would obviously be particularly interesting to investigate 
some more realistic systems using these results. The natural candidate for such an investigation is the <p model with 
short-range interactions in two or more dimensions. At present, only numerical results for the topology of the S w 's 
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are available for this model 10] which do not yet allow to state wheteher the hypotheses of some of the two theorems 
are fulfilled for this model. This is clearly one of the natural lines of development for future investigations. 

We are aware that these sufficient conditions may well be not "optimal" at all, in the sense that it may be possible 
to weaken the hypotheses; moreover, it is possible that in some cases the actual physical relevance of the particular 
phenomenon described here is small, because other mechanisms may be at work |3fJ. Nonetheless, the availability of 
a sufficiency criterion, although restricted to a particular class of transitions and maybe not optimal yet, opens also 
the possibility of using the topological approach to define phase transitions in finite systems, because the topology 
changes in the manifolds T. v which are at the basis of phase transitions do occur also in finite systems: work is in 
progress along this line. Moreover, the phenomenology of "disconnection borders" found in some classical spin systems 
(see e.g. Ref. j31J) might be well related to our present results in the perspective of studying transitional phenomena 
in finite systems. 

In conclusion, we believe that the present work adds some new insight to the topological approach to phase 
transitions, confirming its great potentialities. 
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